Interaction-induced Renormalization of Andreev Reflection 
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We analyze the charge transport between a one-dimensional weakly interacting electron gas and a supercon- 
ductor within the scaling approach in the basis of scattering states. We derive the renormalization group equa- 
tions, which fully account for the intrinsic energy dependence due to Andreev reflection. A strong renormal- 
ization of the corresponding reflection phase is predicted even for a perfectly transparent metal-superconductor 
interface. The interaction-induced suppression of the Andreev conductance is shown to be highly sensitive to the 
normal state resistance, providing a possible explanation of experiments with carbon-nanotube/superconductor 
junctions by Morpurgo et al. [Science 286, 263 (2001)]. 
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The superconducting proximity effect is a well-known phe- 
nomenon, which has motivated a number of theoretical and 
experimental studies since the middle of the last century. The 
low-energy physics of the proximity effect is described by 
Andreev reflection processes 1 1] at the boundary between a 
normal metal (N) and a superconductor (S). In this process 
an electron-like quasiparticle in N is reflected from the NS 
boundary as a hole, thus transferring a double electron charge 
2e into S. The probability of such an event tends to unity in 
the case of an ideal NS interface provided the quasiparticle en- 
ergy e is below the superconducting gap A. Normal reflection 
takes place at non-ideal interfaces due to the Fermi-energy 
mismatch in the superconducting and normal-metal materials 
or due to interface impurities. 

In the elastic theory of electron transport, the NS boundary 
is characterized by energy-dependent quantum-mechanical 
amplitudes: r A (e) for Andreev reflection and r N (e) for nor- 
mal reflection. The differential Andreev conductance of an 
NS junction measured at the voltage bias V is given by |2] 



dl 2e 2 , 



df(e - eV) 
d(eV) 



(l + \r*(er-\r N (er), (1) 



where e is measured with respect to the Fermi energy and 
f{e) is the Fermi distribution function for the temperature T. 
For e ^> A the Andreev amplitude r A vanishes and Eq. ([0 
reduces to the Landauer formula for the conductance. For 
e < A the differential conductance in Eq. Q is solely de- 



termined by r , since \r (e)\ 



i. 



The Landauer formula as well as Eq. Q play a key role 
in the scattering approach, which captures the effects of ge- 
ometry, boundaries, and disorder in nanoscopic samples, but 
ignores inelastic quasiparticle scattering. In this Letter we ap- 
ply Eq. Q to the system depicted in Fig. ^ It consists of 
a superconductor in contact with a one-dimensional electron 
gas with a repulsive interaction, 

(2) 

where U x is a symmetric positive-definite function of x, which 
is assumed to exponentially decay for x exceeding some char- 
acteristic interaction range d. 




x=0 



FIG. 1: The setup consisting of a one-dimensional weakly interact- 
ing electron gas for x < and a superconductor for x > 0. 



It is well-known that electron-electron interactions in a one- 
dimensional Fermi gas strongly modify the spectrum of low- 
energy excitations Q. In this case the elastic amplitudes 

r N,A 
N.A i 



(e) in have to be replaced by the renormalized ones 
r D '^i 6 ) w i tn the cutoff D — max {e, fc^T}. A way to calcu- 
late the renormalized scattering amplitudes was developed in 
Refs. 1 4, 5], where the renormalization group (RG) equation 
for the S-matrix of a single barrier in a one-dimensional inter- 
acting electron gas was derived. The RG results are valid for 
arbitrary barrier transparencies but are restricted to the case 
of weak interactions. Complementary studies based on the 
bosonization technique l6[l7ll8il9 Hl0ll are applicable for arbi- 
trary interaction strengths in the limits of very low or very 
high barrier transparency. The RG method of Refs. fl, y] 
was subsequently extended to account for a resonant energy 
dependence of the bare ^-matrix ifTTl fl2Tl . We demonstrate 
below that the intrinsic energy dependence of Andreev reflec- 
tion, which was disregarded in earlier works Il3lll4ll . is cru- 
cial for understanding the transport properties of interacting 
normal-metal/superconductor nanostructures. 

A possible experimental realization of a quasi-one- 
dimensional electron system is provided by a single-wall car- 
bon nanotube. In recent experiments 1 1511 the differential An- 
dreev conductance was dominated by the interface between 
a carbon nanotube and a superconductor. A strong deviation 
from the predictions of the elastic theory was observed at low 
temperatures and low voltage bias, in which case the Andreev 
conductance was found to be strongly suppressed. The ef- 
fect was shown to disappear with a slight change in the gate 
voltage applied to the nanotube. The major role of the gate 
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voltage in the experimental setup is to govern the nanotube- 
superconductor coupling. Thus, the experiment suggests a 
strong sensitivity of the interaction-induced dip in the An- 
dreev conductance to the normal-state resistance of the NS 
interface. 

Following the RG method of Refs. we deal with 

an effective scattering matrix Sd(e), which acquires an addi- 
tional dependence on the high-energy cutoff D. The renor- 
malization procedure starts at a large D = Do ~ Jwp / d with 
vp the Fermi velocity. The matrix Sd coincides with the bare 
scattering matrix in the absence of interactions 



S Do (e) = S(e) = 



r{e) t'(s) 
t(e) r'(s) 



(3) 



The idea of the S-matrix renormalization is analogous to the 
poor man's scaling proposed by Anderson |16]. One starts 
with the analysis of the first order in U correction to the S- 
matrix, which is regularized by a formal truncation of the 
Fourier series of Hartree and exchange potentials at the large 
momentum cutoff D/hvp. The first order correction is, then, 
proportional to In D je, showing the logarithmic divergence in 
the limit e — ► 0. By differentiating the first order result with 
respect to D one arrives at the RG equation 



dS D (e 



E D - S D (e)T? D S D (e), 



dhxD -wi-d~w» ( 4 ) 

which manifestly conserves the unitarity of the S-matrix fnll . 
The matrix Y>o depends on both Sd{z) and the renormalized 
interaction constants [18] 
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1 - 2ailn (D/D ) 



.92 = a 2 H , (5) 



which makes Eq. (0} a complex non-linear equation. The bare 
values of the interaction constants 



Oil 



dx 



U x e 



2ikpx 



2i:hvF 



Ot 2 



dx ■ 



U, 



2nhvF 



(6) 



quantify backward and forward scattering, correspondingly. 
In order to guarantee the validity of Eq. @ interactions are 
assumed to be weak, i.e. ai,2 <C 1. The RG procedure is 
terminated at D = max{e, ksT}, so that the renormalized 
scattering matrix at zero temperature is given by S e (e). 

For a single resonant tunnel barrier in an interacting quan- 
tum wire the RG equation for the S-matrix can be cast in the 
form of Eq. (0) with MM 



Z D = 



(2g 1 -g 2 )r D (-D) 

(2g' 1 -g' 2 )r' D (-D) 



(7) 



where gi^(D) and g[ 2 {D) refer to the interaction constants 
on the left and the right side of the barrier, correspondingly. 
The matrix Yin takes into account the contribution to the 
Friedel oscillation from electrons with energy —D deep in the 
Fermi see. By lowering the cutoff according to Eq. the 



coherent backward scattering from the electron density oscil- 
lation in the entire energy range is included in the RG proce- 
dure. 

In the presence of superconductivity the dimension of the 
S'-matrix has to be doubled. Therefore, the entries of S{e) in 
the parametrization Q have to be regarded as 2 x 2 matrices 
in Nambu space, for example, 



r(e) 



r N (e) r A {e) 
r A {e) f N {e) 



(8) 



The components of r(e) fulfill an additional electron-hole 
symmetry constraint a y r(—e)*a y = r(e), where a y is a Pauli 
matrix in Nambu space. 

We find that the RG equation for the S'-matrix of the NS 
interface takes the form @ with 



Z D = 



o D 




(9) 



where an is the Nambu matrix 

1 / (2gx - 92)rZ(-D) (gi+g 2 )s 

r N 
' D 

1 



-(r D (-D)+r A ,(D)) 



(2 9l -g 2 )f»(D) " (10) 



(11) 



The lower diagonal block of E/j vanishes due to the absence 
of renormalization on the superconducting side of the inter- 
face. Another important difference from Eq. is the pres- 
ence of the off-diagonal term il lb in <jd, which takes into ac- 
count Friedel oscillations induced by Andreev reflection pro- 
cesses. Even though Andreev reflection can take place for en- 
ergies above the gap, the combination s is non-zero only for 
D < A. Thus, the Andreev renormalization is effective only 
for subgap energies. 

Equation (0J together with Eqs. J9I1 1 i is the main result of 
the present work. The matrix is obtained by differentiating 
the first order perturbation correction to the S-matrix 



d 



dhiD 



dx dy 



Lk, x U X: y(D)Lk, y , (12) 



in the limit et <C D, where Ek = hvpk. The transfer matrix 
Lk :X describes free propagation of electron and hole quasipar- 
ticles 



e i(k F +k)x Q 

g— i(&F— k)x 



(13) 



The first order molecular potential U x . y (D) in Eq. ill 2b de- 
pends on the cutoff and can be decomposed into Hartree and 



exchange terms U x . y 
spondingly, 



H 



u 



jjex 



a z 5{x-y) 



U 



H 



U x x v , that are given by, corre- 



&zU x - z (*l(z)*„(z)), (14) 



~ c^x-y 



2(^(2/)^)} ([^(aO,*^)]) 
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R = 0.001 




FIG. 2: Zero temperature differential Andreev conductance versus 
the voltage bias calculated from Eq. Q with renormalized reflection 
amplitudes. The solid curves result from the numerical solution of 
Eq. I19i with Do = 100A and correspond to different values of the 
normal-state resistance, R = 0.001, 0.01, 0.05, 0.1, 0.2, 0.5, and to 
the particular choice of constants ct\ = 0.1, 012 = 0.2. The dashed 
curves indicate the conductance in the absence of interactions. 



where the square brackets denote the commutator. The op- 
erators ^ in Eq. (I14> have to be taken in the scattering basis 



dk 



(V( fc 



F+k)x 



-k)x^ 



(15) 



where the momentum integration is restricted to the interval 
(—D/hv-p, D/hvp), leading to the cutoff dependence of U x ,y 
in Eq. The operators a a % and b a k are related by the 

reflection matrix JHl 



?Tfc 



& 1- 



(16) 



Vfc L V/c' 



The angular brackets in Eq. dl4t are defined by 
f(£k)5aa'fikk>- Calculating ap from Eq. fl!2i in the limit 
et,T«Dwe obtain the result l fl"0l . 

Both the Hartree and exchange potentials in Eq. dl4i are 
hermitian, x = C^r.j/i and fulfill an electron-hole symme- 
try, J7J x = U x ,y> for any £>. Thus, both the unitarity as well 
as the electron-hole symmetry of the S'-matrix are conserved 
under the RG flow (g). 

The bare S'-matrix for Andreev reflection has a specific en- 
ergy dependence on the scale A. For an ideal NS interface the 
bare Andreev reflection amplitude is given by 



= e/A - V(e/A) 2 - 1, e > 0, 



(17) 



where the positive branch of the square root has to be chosen 
for e < A and the relation 7_ e = —7* extends the definition 
to negative energies. Furthermore, the reflection matrix l|8} 
of a non-ideal interface can be conveniently parameterized as 
d 



'Mi 



(1 ~ 7e 2 )Vfi 

1-HR ' 



(18) 



where we assume the phases 6x,% and the normal-state reflec- 
tion coefficient Re (0, 1) to be energy independent. 

Due to the simple structure of in Eq. it suffices to 
consider the reflection block of Eq. © 



dr D (e) 
dhxD 



rD{£)^ D r D (e) 



(19) 



In what follows we focus on the physically relevant limit of 
short-range interactions d < £ or, equivalently, Do > A, 
where £ = Hvf / A is the superconducting coherence length. 
We start from the analysis of Eq. (1191 for the case of an ideal 
interface, R — 0, and parameterize f^(e) = exp (— i4>£>{e)), 
where the phase 4>d{£) is real for e < A. The Andreev con- 
ductance below the gap is therefore not affected by interac- 
tions, while the Andreev phase 4>d(s) strongly deviates from 
its elastic value arccose/A. Such an Andreev phase renor- 
malization is absent in Refs. QUI. 
With the substitution 



^d(e) = arccos ■ 



r/A 



(20) 



cosh it £> — dl — (e/A) 2 sinhw£> 
Eq. d!9t is reduced to the differential equation for a function 



of one variable 

duo 



(gi(D) + g 2 {D)) tanh(u D - w D ), 



(21) 



with the initial condition ua = 0, and with the notation wd — 
arctanh \/l — (D/A) 2 . The solution to Eq. (12 II in the limit 
D <C A can be approximated by 

ujj = (a.2— oii/2) In ^ + m f 1 + a m J^j + cons t-, (22) 

where the parameter a = 2ai / ( 1 + 2a x In D / A) has a log- 
arithmic dependence on the initial cutoff Do ~ hv F / d. From 
Eq. (I20> we obtain 



cos<t) D {e) -^(jj 



12 —Ot\j% 



1 



3/4 



(23) 



where = max{e, fcgT} < A. At zero temperature we let 
D = e in Eq. d23i and interpret the result as an additional 
suppression of the electron-hole coherence at finite e. The ef- 
fect of the Andreev phase renormalization can be seen most 
explicitly in the supercurrent or in measurements of the den- 
sity of states. A detailed analysis of these quantities, though, 
is beyond the scope of the present work. We shall stress, how- 
ever, that the renormalization d23l > originates in the energy de- 
pendence of the bare scattering amplitudes, which has been 
ignored in the previous studies of the Josephson effect and of 

the density of states ilia. 

For a non-ideal interface, R 7^ 0, Eq. dl9t reduces to a set 
of four coupled RG equations for the variables 7, R and B\ 2 
in the parameterization i ll 8b . It is remarkable that for energies 
above the gap only one equation remains 



dR D 
dlxiD 



{2 9l - g 2 )R D {l - R D ) 



1 



1 - 1IR1 



(24) 
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with the initial condition Ro = R and with the function 7^ 
given by the bare Andreev amplitude d!7i . Both the ampli- 
tude 7 £ and the phases 6\ t 2 are not renormalized and do not 
acquire any D dependence as far as D > A. It, then, follows 
from Eq. (II 81 that r N (A) = 0, hence the differential Andreev 
conductance at eV = A equals its elastic value 4e 2 /h. On 
the other hand, renormalization (I24> reduces to the result of 
Refs. 1 4, 5] in the limit A — * and yields the well-known con- 
ductance suppression at low energies. Thus the interactions 
tend to sharpen the non-monotonic behavior of the Andreev 
conductance near eV = A in accordance with Ref. 12(111 . 

For energies below the gap, Eq. ( 1191 is equivalent to a 
joint renormalization of the three variables cj>, R and 9% in 
the parameterization (II 81 . where is the Andreev phase, 
7 = exp (— i(j>). Thus, below the gap the renormalization 
cannot be reduced to that of the single parameter R in con- 
tradiction to the phenomenological description of Ref. 12 ill . 
The phase O2 is not renormalized in the entire energy range, 
which is related to the fact that it can be removed by an appro- 
priate gauge transformation. The Andreev conductance below 
the gap depends only on the absolute value of the amplitude 
r^3 (e). One can prove from the explicit form of the RG equa- 
tions that the energy dependence of |rg(e)|, unlike that of 
r^(e), can be disregarded in the limit e <C A. Therefore, we 
obtain from Eq. (d ESQ 



A\2 



dlnD 



= 2{2a 2 -ai)\r£\ 2 (l-\r*\ 2 ), D < A, (25) 



which depends only on the invariant combination 2a 2 — ol\ of 
constants, which is positive for any repulsive interaction. The 
general solution of Eq. J25l > can be written as 



4? 



l + b(A/D) 



2(2a 2 -Qi) ' 



(26) 



where the coefficient b has to be determined from the solution 
of Eq. dl9l for e, D ~ A. For a nearly ideal interface, R -C 1, 
we find b oc R, which allows to estimate the width of the 
zero bias anomaly in the differential Andreev conductance as 
e c ~ Ai? 1 ^ 402-2 " 1 ). 

We illustrate our findings in Fig. [2] The differential An- 
dreev conductance at zero temperature is calculated from 
Eq. Q with the renormalized amplitudes r^' A (e), which are 
found by the numerical solution of Eq. il9\ . The interaction 
induced dip in the Andreev conductance has a width eV c ~ e c 
and is strongly sensitive to the normal-state resistance of the 
interface. Such sensitivity was indeed observed in the experi- 
ments by Morpurgo et al. flUl . where the normal-state trans- 
parency of the NS junction was controlled by means of a gate 
voltage applied to the nanotube. A small change in the gate 
voltage gave rise to a slight improvement of the normal-state 
transparency of the junction, which had a drastic effect on the 
Andreev conductance. 



In conclusion we have derived and analyzed the energy- 
dependent RG equations for the scattering matrix of an inter- 
acting normal-metal/superconductor interface. Our approach 
takes into account the intrinsic energy dependence of scatter- 
ing at the NS interface and is readily generalized to SNS struc- 
tures and to the case of a quasi-one-dimensional interacting 
normal metal. The effects of interaction on the Andreev con- 
ductance, but not on the Andreev phase, are shown to vanish 
in ideal NS junctions. Our results qualitatively explain recent 
experiments fl5ll with carbon nanotubes. 
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